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A method of and apparatus for processing Seismic data 

The present invention relates to a method of processing multi-component seismic data. 
It particularly relates to a method of processing seismic data to determine a calibration 
filter that calibrates one component of the seismic data relative to another component of 
the seismic data. The invention further relates to an apparatus for processing seismic 
data. 

Figure 1 is a schematic view of a seismic surveying arrangement. In this figure the 
surveying arrangement is a marine surveying arrangement in which seismic energy is 
emitted by a seismic source 1 that is suspended within a water colxmm 2 firom a towing 
vessel 3. When the seismic source 1 is actuated seismic energy is emitted downwards . 
and is detected by an array of seismic receivers 4 disposed on the seafloor 5. (As used 
herein the term "seabed" denotes the earth's interior, and the term "seafloor" denotes 
the surface of the seabed.) 

Many seismic surveys now use multi-component receivers that record two or more 
components of the seismic energy incident on the receiver. For example a 3-component 
(3-C) seismic receiver contains three orthogonal geophones and so can record the x-, y- 
and z-components of the particle motion at the receiver (the particle motion may be the 
particle displacement, particle velocity or particle acceleration or even, in principle, a 
higher derivative of the particle displacement). In a marine seismic survey a 4- 
component (4-Q seismic receiver can alternatively be used. A 4-C receiver contains a 
pressure sensor such as a hydrophone in addition to three orthogonal geophones and so 
can record the pressure of the water column (which is a scalar quantity) in addition to 
the X-, y- and z-components of the particle motion. 



Many different paths exist by which seismic energy may travel fi"om the source 1 to a 
receiver 4 in the seismic surveying arrangement of Figure 1. A number of paths are 
indicated schematically in Figure 1, 
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The path 6 shown in Figure 1 is known as the "direct path". Seismic energy that travels 
along the direct path 6 travels from the source 1 to a receiver 4 essentially in a straight 
line without undergoing reflection at any interface. 

Path 7 in Figure 1 is an example of a "water layer multiple path". Seismic energy that 
follows a water layer multiple path propagates wholly within the water column 2, but 
undergoes one or more reflections at the surface of the water column and/or the seafloor 
5 so that the seismic energy passes through the water column more than once. The 
water layer multiple path 7 shown in Figure 1 involves one reflection at the seafloor 5 
and one reflection at the surface of the water column, but many other water layer 
multiple paths exist. 

The path 8 in Figure 1 is an example of a "critical refraction path". Seismic energy that 
follows the path 8 propagates dov^wards to the seafloor 5, and penetrates into the 
earth's interior 10 (ie into the seabed). The seismic energy continues propagating 
downwardly, until it reaches a boundary 1 1 between two layers of the earth that have 
different acoustic impedance. The seismic energy undergoes critical refraction, 
propagates along the boundary 1 1, before eventually being refracted upwards towards 
the receiver 4. Critical refraction may also occur at the water-seabed interface, and 
downwardly propagating seismic energy that is refracted in this way will propagate 
along the water-seabed interface and will then propagate upwardly into the water 
column. 

The path 9 shown in Figure 1 is known as a "primary reflection path". Seismic energy 
that follows the primary reflection path 9 propagates downwards through the water 
column, is refracted at the seafloor 5, and propagates downwardly through the earth's 
interior. The seismic energy is refracted at the boundary 1 1, but is not critically 
refracted and so continues to propagate downwardly into the earth. It eventually 
undergoes reflection at a geological structure 12 that acts as a partial reflector of seismic 
energy, and the reflected seismic energy is, after fiirther refraction as it passes upwardly 
through the boundary 1 1, incident on the receiver 4. The general intent of a seismic 
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survey is to make use of the seismic energy that follows the primary reflection path in 
order to obtain information about the interior structure of the earth. 

Seismic energy acquired at a receiver may contain upwardly and/or downwardly 
propagating seismic energy depending on the location of the receiver and on the event. 
For example seismic energy that travels along the critical refraction path 8 shown in 
Figure 1 will, when it is incident (travelling upwardly) on the water-seabed interface, be 
partly transmitted into the water column and partially reflected back into the seabed 10. 
Thus, a critical refraction event will consist purely of upwardly propagating seismic 
energy above the seafloor 5, but will contain both upwardly and downwardly 
propagating seismic energy below the seafloor 5. As another example, seismic energy 
that travels along the direct path 6 shown in Figure 1 will, when incident on the water- 
seabed interface 5, be partially transmitted into the seabed and partially reflected back 
into the water column. Hence, the direct event will contain both upwardly and 
downwardly propagating seismic energy above the seafloor, but will contain only 
downwardly propagating seismic energy below the seafloor. It is therefore often of 
interest to decompose the seismic data acquired at the receiver 4 into an up-going 
constituent and a down-going constituent, above or below the seafloor 5, For example, 
in a 4-C seismic survey it may be of interest to decompose the pressure and the vertical 
particle velocity recorded at the receiver into their up-going and down-going 
constituents above the seafloor. Various filters that enable decomposition of seismic 
data into up-going and down-going constituents have been proposed. One example can 
be found in K.M. Schalkwijk et al, "Application of Two-Step Decomposition to Multi- 
Component Ocean-Bottom Data: Theory and Case Study", J. Seism. Expl. Vol. 8 
pp261-278 (1999), and states that the down-going and up-going constituents of the 
pressure just above the seafloor may be expressed as follows: 
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where P is the pressure acquired at the receiver, F is the up-going constituent of the 
pressure above the seafloor, is the down-going constituent of the pressure above the 
seafloor,/is the frequency, k is the horizontal wavenumber, Z is the vertical particle 
velocity component acquired at the receiver, ji? is the density of the water, and q is the 
vertical slowness in the water layer. 

As can be seen, the expressions in equation (1) require two of the components of 
seismic data recorded at the receiver to be combined. These filters are an example 
where it is necessary to combine two components of the acquired seismic data. It may 
also be necessary to combine two or more components of the. acquired seismic data in 
order to decompose the acquired seismic data into p-wave and s-wave (pressure-wave 
and shear-wave) components, or to remove water level multiple events from the seismic 
data. 

One problem in combining different components of the seismic data acquired at a 
receiver is that the different components of the seismic data may not be correctly 
calibrated against one another. This is particularly the case where the two components 
that are being combined are, as in equation (1), the pressure and the vertical particle 
velocity. There are usually differences in coupling or impulse response between the 
hydrophone used to acquire the pressure and the geophone used to acquire the vertical 
particle velocity. It is necessary to calibrate the data for these differences before the 
pressure and vertical particle velocity can be combined. This may be done by 
developing a calibration filter that compensates for the differences in coupling and 
impulse response between the hydrophonie and the vertical geophone. 

Schalkwijk et al, and others, have suggested that the calibration problem can be 
addressed by assuming that one component of the seismic data has been correctly 
recorded, and calibrating the other component of the seismic data against the component 
that is assumed to be correctly recorded. In general, it is assumed that the hydrophone 
is well coupled, so that the pressure recording is taken to be correct. The vertical 
component of the particle velocity is then calibrated against the pressure to compensate 
for coupling and impulse response differences between the hydrophone and the vertical 
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geophone. Schalkwijk et al therefore proposed that equation (1) above should be 
modified by applying a calibration filter to the vertical particle velocity. They proposed 
that the equation given above for the down-going constituent of the pressure above the 
seafioor should be modified to read as follows: 

rif,k)^\pif,k)^aif)^^Zif,k). _ (2) 

In equation (2) a(f) represents a fi-equency-dependent calibration filter. The remaining 
terms in equation (2) have the same meaning as in equation (1). 

The method proposed by Schalkwijk et al for detennining the calibration filter a(f) is to 
minimise the energy of the down-going pressure constituent above the seafloor for a 
portion of the seismic data that contains only primary reflections. Seismic energy 
travelling along a primary reflection path is propagating upwardly just above the 
seafloor at the receiver position, so that the down-going constituent of the pressure just 
above the seafloor should be zero for data that contains only primary reflections. 
Schalkwijk proposed that the calibration filter that minimises the energy of the down- 
going pressure in a window containing only primary reflection events can be found 
using a least squares method. Once the calibration filter a(f) has been determined in this 
way, it is applied to the entire data set for the vertical particle velocity. 

The existence of various paths of seismic energy from the source to the receiver means 
that the data acquired at the receiver in a real seismic survey will contain events 
corresponding to more than one possible path. These events will occur at different 
times after the actuation of the seismic source 1, as different paths of seismic energy 
have different associated travel times. Figure 2 is a schematic illustration of seismic 
data that might be acquired at the receiver 4, and it shows the amplitude of seismic 
energy recorded at the receiver 4 as a fimction of the time since the actuation of the 
source 1 . Figure 2 illustrates a direct event 13, coixesponding to the direct path 6, a 
critical refraction event 14 corresponding to the critical refi-action path 8, a primary 
event 15 corresponding to the primary reflection path 9, and a water layer multiple 
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event 16 coTresponding to the water layer multiple path 7. (In practice, data acquired at 
a receiver will contain a plurality of primary reflection events fix)m different geological 
structures, a plurality of critical refraction events, and a plurality of water level mxiltiple 
events arising from different water level multiple paths. Only one event of each type is 
shown in Figure 2 for simplicity of explanation.) In order to apply the method of 
Schalkwijk et al to determine the calibration filter, data in a time window that contains 
only the primary event 15, such as the time window A shown in Figure 2, must be 
selected. 

The present invention provides a method of processing multi-component seismic data 
obtained from seismic signals propagating in a medium, the method comprising the 
steps of: selecting a first portion of the seismic data containing only events arising from 
critical refraction of seismic energy; and determining a first calibration filter from the 
first portion of the seismic data, the first calibration filter being to calibrate a first 
component of the seismic data relative to a second component of the seismic data. 

The method proposed by Schalkwijk et al has the disadvantage that the time window 
containing only primary reflection events has to be picked manually. The primary 
reflection events are not the first events acquired at the receiver following actuation of 
the source, and so cannot be picked automatically. A further disadvantage is that in 
some cases, for example if the seismic source has a long signature, it may be hard to 
distinguish between the direct arrival and the primary reflection events, so that it may be 
difficult to isolate the correct events. The direct event contains downwardly 
propagating seismic energy so that use of a time window that inadvertently included the 
direct event would not give correct results for the calibration filter, since the method for 
determining the calibration filter assumes that the selected data contains only up-going 
energy. A further problem with the method of Schalkwijk et al is that in shallow water 
the water layer multiple events may arrive at substantially the same time as the primary 
reflection events, and this again makes it difficult to pick a time window that includes 
only the primary reflection events. 
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The present invention makes use of the fact that the critical refraction events consist 
only of up-going seismic energy just above the seafloor. Thus, selecting a time window 
that contains only one or more critical refraction events makes it possible to determine 
the calibration filter a(j) by the technique of minimising the energy of the down-going 
pressure just above the seafloor in that time window. 

The method of the invention is particularly advantageous when applied to long offset 
data. As is shown in Figure 3, as the offset (that is, the horizontal distance between the 
source and the receiver) increases, the arrival time of the first critical refraction event 
increases more slowly than does the arrival time of the direct event. For offsets greater 
than Oi the first arrival at the receiver is not the direct event, but is the critical refraction 
event. That is, at long offsets the critical refraction event 14 in Figure 2 (two critical 
refraction events are shown in Figure 3) will arrive before the direct event 13 and will 
be the first arrival at the receiver When the invention is applied to data having a source- 
receiver offset sufficiently large for the first event acquired at the receiver to be a 
critical refraction event, it is possible to use a time window that covers only the first 
event acquired at the receiver - and this makes it possible to use an automatic picking 
method to detemiine the time window. If several critical refraction events arrive at the 
receiver before the direct event, as for far offsets in Figure 3 where two critical 
refraction events anive before the direct event arrives, then all these critical refraction 
events may be included in the time window. 

Figure 4 is a schematic illustration that corresponds to Figure 2, but illustrates the 
arrival times of the events at an offset that is sufficientiy large such that the first arrival 
is a critical reflection event. In this case, the invention may be applied by selecting the 
time window B which includes the critical refraction event only, and minimising the 
energy of the down-going pressure above the seafloor in this time window. 

A fiirther advantage of the invention is that the method may be applied to seismic data 
acquired in shallow waters. Although water layer multiple events in seismic data 
acquired in shallow water may coincide with primary reflection events, they do not 
coincide with critical refraction events. Choosing a time window that includes only the 
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critical refraction event therefore ensures that the time window cannot contain water 
layer multiple events. The invention also overcomes the problems that arise when a 
seismic source having a long source signature is used. 

A preferred embodiment of the invention comprises the further steps of selecting a 
second portion of the seismic data containing only events arising from primary 
reflection of seismic energy and determining a second calibration filter from the second 
portion of the seismic data, the second calibration filter being to calibrate the first 
component of the seismic data relative to the second component of the seismic data. It 
may comprise the further step of determining a wavenumber-dependent calibration filter 
from the first calibration filter and the second calibration filter. 

A further problem related to the method proposed by Schalkwijk et al is that the correct 
calibration filter a(f) may well be dependent on the wavenimiber as well as on the 
frequency. The calibration filter proposed by Schalkwijk, however, is dependent only 
on frequency and, furthermore, is derived purely from seismic data at low 
wavenumbers. In an embodiment of the present invention, the filter obtained from the 
critical refraction events is combined with a filter obtained from primary reflection 
events, and a wavenumber-dependent filter is obtained from the two individual filters. 
The wavenumber-dependent filter may be obtained by, for example, inteipolation 
between the filter derived from the critical refraction events and the filter derived from 
the primary reflection events. 

A second aspect of the present invention provides a method of processing multi* 
component seismic data obtained from seismic signals propagating in a medium, the 
method comprising the step of selecting a first portion of the seismic data corresponding 
to a first wavenumber range; determining a first calibration filter from the first portion 
of the seismic data; selecting a second portion of the seismic data corresponding to a 
second wavenumber range different from the first wavemmfiber range; determining a 
second calibration filter from the second portion of the seismic data; and determining a 
wavenumber-dependent calibration filter from the first calibration filter and the second 
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calibration filter, tiie wavenumber-dq)endent calibration filter being to calibrate a first 
component of the seismic data relative to a second component of the seismic data. 

A third aspect of the invention provides a method of processing multi-component 
seismic data obtained from seismic signals propagating in a medium, the method 
comprising the steps of: selecting a first portion of the seismic data in which the first 
" arrival contains only upwardly propagating seismic energy above the seafloor; and 
determining a first calibration filter from the first portion of the seismic data, the first 
calibration filter being to calibrate a first component of the seismic data relative to a 
second component of the seismic data. 

The invention may be applied to any event that is the first arrival and that contains only 
upgoing energy above the seafloor. For example, at far offsets the first arrival may be 
an event which is not a critical refraction event but which nevertheless contains only 
upgoing energy above the seafloor - such as, for example, a wave that was trapped in a 
thin subsurface layer of the seabed - and the invention may be applied to such events. 

The invention may further comprise the step of calibrating the first component of the 
seismic data using the first calibration filter or using the wavenumber-dependent 
calibration filter. 

A fourth aspect of the present invention provides a method of seismic surveying 
comprising the steps of: actuating a source of seismic energy; acquiring seismic data at 
a receiver spatially separated from the source; and processing the seismic data by a 
method as defined above. 

A fifth aspect of the present invention provides an apparatus for processing multi- 
component seismic data to determine a calibration filter to calibrate a first component of 
the seismic data relative to a second component of the seismic data, the apparatus 
comprising: means for selecting a first portion of the seismic data containing only 
events arising from critical refraction of seismic energy; and means for determining a 
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first calibration filter from the first portion of the seismic data. The apparatus may 
comprise a programmable data processor. 

A sixth aspect of the invention provides an apparatus for processing multi-component 
seismic data to determine a calibration filter to calibrate a first component of the seismic 
data relative to a second component of the seismic data, the apparatus comprising: 
means for selecting a fu-st portion of the seismic data in which the first arrival contains 
only upwardly propagating seismic energy above the seafloor; and means for 
detemiining a first calibration filter firom the first portion of the seismic data. 

A seventh aspect of the invention provides an apparatus for processing multi-component 
seismic data to determine a calibration filter to calibrate a first component of the seismic 
data relative to a second component of the seismic data, the apparatus comprising: 
means for selecting a first portion of the seismic data corresponding to a first 
wavenumber range; means for detemiining a first calibration filter from the first portion 
of the seismic data; means for selecting a second portion of the seismic data 
corresponding to a second wavenumber range different from the first wavenumber 
range; means for determining a second calibration filter from the second portion of the 
seismic data; and means for determining a wavenumber-dependent calibration filter 
from the first calibration filter and the second calibration filter. 

The apparatus may comprise a programmable data processor. 

An eighth aspect of the present invention provides a storage medium containing a 
program for an apparatus as defined above. 

The invention also provides a method of determining a first calibration filter for 
calibrating a fu:st component of multi-component seismic data relative to a second 
component of the seismic data, the method comprising the steps of: selecting a first 
portion of the seismic data containing only events arising from critical refraction of 
seismic energy; and determining the first calibration filter from the first portion of the 
seismic data. 
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The invention also provides a method of determining a wavenumber-dependent 
calibration filter for calibrating a first component of multi-component seismic data 
relative to a second component of the seismic data, the method comprising the steps of: 
selecting a first portion of the seismic data corresponding to a first wavenumber range; 
determining a first calibration filter from the first portion of the seismic data; selecting a 
second portion of the seismic data corresponding to a second wavemmaber range 
different fi'om the first wavenumber range; determining a second calibration filter from 
the second portion of the seismic data; and determining a wavenumber-dependent 
calibration filter from the first calibration filter and the second calibration filter. 

The invention also provides a method of determining a first calibration filter for 
calibrating a first component of multi-component seismic data relative to a second 
component of the seismic data, the method comprising the steps of: selecting a first 
portion of the seismic data in which the first arrival contains only upwardly propagating 
seismic energy above the seafloor; and determining a first calibration filter from the first 
portion of the seismic data. 

Preferred embodiments of the present invention will now be described by way of 
illustrative example with reference to the accompanying figures in which: 

Figure 1 is a schematic illustration of a seismic survey; 

Figure 2 is a schematic illustration of the seismic energy acquired at a receiver in the 
seismic survey of Figure 1; 

Figure 3 is a schematic illustration of the variation of arrival time of seismic energy as a 
function of offset between the source and the receiver; 

Figure 4 is a schematic illustration of the seismic energy acquired in the seismic survey 
of Figure 1 at long offsets, illustrating a method of the present invention; 
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Figure 5 is a schematic illustration of pressure recorded at a receiver in the seismic 
surveying arrangement of Figure 1; 

Figures 6 and 7 illustrate the up-going and down-going constituents of pressure above 
the seafloor obtained from the pressure data shown in Figure 5 according to a prior art 
approach; 

Figures 8 and 9 illustrate the up-going and down-going pressure constituents above the 
seafloor obtained from the pressure data of Figure 5 according to a method of the 
present invention; 

Figure 10 is a schematic block flow diagram of a method of the present invention; and 

Figure 1 1 is a block schematic diagram of an apparatus according to the present 
invention. 

Figure 5 illustrates typical pressure data recorded at a 4-C receiver in a seismic survey 
such as the survey shown in Figure 1 . The x-axis in Figure 5 indicates the offset 
between the source and the receiver, and the y-axis indicates the time after actuation of 
the seismic source. The data are common receiver data and were acquired using a 
single receiver and a linear array of sources deployed with a spacing of 50m between 
each pair of adjacent sources. Each trace represents the pressure acquired at the 
receiver when one source is actuated, with the amplitude of the acquired pressure being 
in the x-direction. 

It should be noted that different receivers in an array may well have different coupling, 
different instruments responses etc, even if all the receivers are nominally identical to 
one another. The calibration fiher required for data acquired at one receiver in a 
receiver array is therefore likely to be different from the calibration filter required for 
data acquired at another receiver in the array. The invention is therefore preferably 
applied to common receiver gathers and a separate calibration filter is determined for 
each common receiver gather. 
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The pressure data shown in Figure 5 contains a large number of seisnaic events. The 
event labelled 13 is the direct wave, and it will be seen that this is the first arrival for 
offsets having a magnitude of up to approximately 1,000m. The event labelled 14 is a 
critical refiraction event, and it will be seen that this is the first arrival for offsets having 
a magnitude significantly greater than 1,000m. 

Figures 6 and 7 illustrate the up-going constituent above the seafloor (Figure 6) and the 
down-going constituent above the seafloor (Figure 7) of the pressure shown in Figure 5 
obtained using the filters given in equation (1) above. That is, the up-going and down- 
going constituent shown in Figures 6 and 7 were obtained on the assximption that the 
pressure data and the vertical particle velocity data (not illustrated) were correctly 
calibrated to one another. Inspection of Figures 6 and 7 shows that this assumption is 
incorrect. In particular, the critical refraction event 14 contains only up-going energy 
above the seafloor and so should appear only in the up-going pressure constituent and 
should not appear in the down-going pressure constituent. It will, however, be seen that 
the up-going critical refraction event has leaked through into the down-going pressure 
constituent shown in Figure 7, and this indicates that the calibration is unsatisfactory. 

According to the present invention, a calibration filter is determined from the critical 
refraction event 14. As noted above, for traces acquired at a source-receiver offset 
having a magnitude well above 1000m, the critical refraction event is the first event 
acquired at the receiver, and is well-separated from the subsequent event. It is therefore 
possible for such traces to define a time-offset window that includes only the first 
critical refraction event, and so includes only up-going energy. 

One suitable time-offset window of data is illustrated in Figure 5, as the region C. It 
will be seen that this region includes traces acquired at an offset of between -3000m to 
approximately -2100m. For each trace in this offset range the region C defines a time 
window that includes only the first refraction event (which is the first arrival in each of 
the selected traces). It will be noted that the centre point of the time window for a 
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particular trace is not constant between traces but increases with increasing magnitude 
of offset. 

The calibration filter for the vertical velocity component is then calculated on the 
assumption that the energy in the selected portion C of the data should contain only up- 
going energy. The calibration filter may be determined m any suitable way. In 
particular, the calibration fiUer a(f) may be determined by finding the calibration filter 
that minimises the energy of the down-going pressure constituent using a least squares 
process, as in the method of Schalkwijk et al. Once the appropriate cahbration filter a(j) 
has been determined, revised filters for determining the up-going and down-going 
constituents of the pressure above the seafloor can be determined using equation (2), or 
in general the filter a(f) may be applied to the entire gather of the vertical component 
data, and the calibrated vertical component data can then be used as an input to any 
process requiring a combination of the vertical component with any other seismic 
components. 

Another suitable portion of data exists in the corresponding region for positive offsets in 
the range 2100 to 3000m. One possible implementation of the method would be to use 
both these regions, by defining a second region, analogous to the region C in Figure 5, 
for offsets in the range +2 100m to +3000m and determining a second calibration filter. 
The two filters determined fi^om the two windows may then be averaged. This will 
however not be possible for all data sets since a receiver gather does not necessarily 
have the same amount of positive and negative offsets, and hence a region with a clearly 
sq)arated critical refiraction event may be present only for either positive or negative 
offsets. 

Figures 8 and 9 illustrate the results of decomposing the pressure data of Figure 5 into 
its up-going and dovm-going constituents above the seafloor using filters of the type 
given in equation (2) above, and with the calibration filter a(f) determined fi-om the 
seismic data in the region C of Figure 5. It will be noted that the critical refi-action 
event 14 appears predominantly in the up-going pressure constituent of Figure 8, and is 
almost completely absent fi-om the down-going pressure constituent of Figure 9. This 
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illustrates that the decomposition of Figures 8 and 9 is significantly more accurate than 
the decomposition of Figures 6 and 7, since the critical refiraction event is expected to 
occur only in the up-going pressure constituent. 

It will also be noted that the primary reflection event is stronger in the up-going 
pressure constituent of Figure 8 than in the up-going pressure constituent of Figure 6. 
This suggests that the calibration filter found from the critical refraction event at long 
offsets is also applicable at low offsets. 

Figure 10 is a block flow diagram illustrating one embodiment of the method of the 
present invention. 

Initially, at step 17, seismic data is acquired. This may be, for example, acquired in a 
survey of the type shown in Figure 1 . 

The invention may altematively be applied to pre-existing seismic data. Step 17. may 
therefore be replaced by the alternative step 18 of retrieving pre-existing seismic data 
from storage. 

At step 19 a suitable offset range is selected. In the example described above with 
reference to Figure 5, step 19 consists of selecting the offset range from -3,000m to - 
2100m. 

At step 20, the first arrival of seismic energy for each trace in the selected offset range is 
determined (this may be thought of a selecting a time window for each trace, and so 
defining an offset-time window). Assuming that the offset range has been selected 
correctly in step 19, the first arrival in each trace in the selected offset range will be a 
critical refiraction event such as the event 14. Since the event is the first event in each 
trace, step 20 maybe carried out using an automatic picking method, although it may 
altematively be performed manually. 
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At step 21, a calibration filter is determined that is the best fit to the data in the selected 
offset range and time window. This is done by calculating the down-going pressure 
constituent above the seafloor from the pressure and vertical particle velocity recorded 
at the receiver using equation (2), and finding the calibration filter that minimises the 
energy in the down-going constituent of the pressure. 

At step 22 the filter a(f) is applied to all the desired traces of the vertical component of 
the seismic data acquired at step 17 or retrieved fi-om storage at step 18. 

At step 23 the calibrated vertical component data is then used as input into any process 
requiring a combination of several seismic components. For example, filters for 
detemiining the up-going and down-going constituents of the pressure above the 
seafloor may be determined, using equation (2) and the calibration filter determined at 
step 21. 

If desired, steps 22 and 23 may be omitted. In this case the calibration fiher determined 
at step 21 may be output for display or stored for subsequent use. 

In an alternative embodunent of the invention, a wavenumber-dependent filter is 
determined by combining the approach of Schalkwijk et al with the present invention. 
In this embodiment, a calibration filter is determined firom critical refi-action events 
occurring at long offsets, as described above with reference to steps 17 to 21 of Figure 
10. A second calibration filter is then determined firom traces in which it is possible to 
defme a time-offset window that contains only primary reflection events. A calibration 
filter is determined for these traces in the manner described by Schalkwijk et al. A 
suitable region of data for obtaining this filter is indicated on Figure 5 as region D. 

The calibration filter determined from critical refraction events at long offset and the 
calibration filter determined from primary reflection events at low offset are then 
combined to produce a wavenumber-dependent calibration filter. The filters may be 
combined using an interpolation technique to determine the filter to be applied at a 
given offset. 
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In this embodiment, step 22 of Figure 10 is replaced by the step of calibrating the 
vertical component using the wavenumber-dependent calibration filter. Alternatively 
steps 22 and 23 may be omitted, and the wavenumber-dependent calibration filter can 
be output or stored for fiiture use. 

An alternative way to obtain a wavenumber-dependent calibration filter is to compute a 
calibration filter for each separate trace in the offset range selected at step 19. In this 
alternative embodiment steps 20 and 21 are performed on each trace (or on a plurality of 
selected traces) in the offset range selected at step 19 so that calibration filters are 
determined for several different wavenumbers. Alternatively, the traces in the offset 
range selected at step 19 can be grouped, and a calibration filter can be determined for 
each groiq) of traces, for example using a least squares method. Again, this results in 
calibration filters for several different wavenumbers. 

Once calibration filters have been obtained for several different wavenumbers, it is 
possible to interpolate between and/or extrapolate fi-om these calibration filters to obtain 
a wavenumber-dependent calibration filter. This method would, however, only work 
well for a time-offset window containing only primary reflections (i.e., the window D in 
Figure 5), since the vertical slowness is constant for the refracted event. The 
wavenumber-dependent calibration filter may be used immediately, or may be output or 
stored for fiiture use. 

A fiirther altemative method is to define time-offset windows around several refraction 
events of different vertical slownesses, and determine a plurality of calibration filters 
(one calibration filter can be obtained firom data in each window). A wavenumber- 
dependent calibration filter can be obtained by interpolation between and/or 
extrapolation fi-om these calibration filters. If desired, one or more calibration filters 
determined fi-om a time-offset window containing only primary reflections can also be 
used in the interpolation and/or extrapolation. The wavenumber-dependent cahbration 
filter again may be used immediately, or may be output or stored for future use. 
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It will be noted in Figure 5 that the primary reflection event is obscured by other events 
at long offsets. It will therefore be extremely difficult to compute a reliable calibration 
filter at long offsets using the method of Schalkwijk et al, owing to the difficulty of 
determining a time window that contains only primary reflection events. Furthermore, 
even if a time window that contained only primary reflection events could be 
determined for the long offset traces in Figure 5, this could only be done by a manual 
picking method and could not be automated. 

The invention has been described above with reference to a calibration filter that 
calibrates the vertical particle motion with regard to the pressure, on the assumption that 
the pressure has been accurately recorded. The invention is not limited to this, however, 
and in principle could be used to determine a calibration filter that calibrates the 
pressure with regard to the vertical particle motion, on the assumption that the vertical 
particle motion has been accurately recorded. 

Figure 1 1 is a schematic block diagram of an apparatus 34 according to the present 
invention. The apparatus is able to carry out a method according to the present 
invention. 

The apparatus 34 comprises a progranmiable data processor 27 with a program memory 
28, for instance in the form of a read only memory ROM, storing a program for 
controlling the data processor 27 to process seismic data by a method of the invention. 
The apparatus further comprises non-volatile read/write memory 29 for storing, for 
example, any data which must be retained in the absence of power supply. A "working" 
or "scratchpad" memory for the data processor is provided by a random access memory 
(RAM) 30. An input device 31 is provided, for instance for receiving user commands 
and data. An output device 32 is provided, for instance for displaying information 
relating to the progress and result of the method. The output device may be, for 
example, a printer, a visual display unit or an output memory. 

Seismic data for processing may be supplied via the input device 31 or may optionally 
be provided by a machine-readable store 33. 
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The program for operating the apparatus and for performing a method as described 
hereinbefore is stored in the program memory 28, which may be embodied as a semi- 
conductor memory, for instance of the well-known ROM type. However, the program 
may be stored in any other suitable storage medium, such as magnetic data carrier 28a 
(such as a "floppy disc") or CD-ROM 28b. 



